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ABSTRACT 

We study the influence of jet structure on the curvature effect in a GRB-lightcurve. Using a simple 
model of jet emissivity, we numerically calculate lightcurves for a short flash from a relativistic outflow 
having various profiles of the Lorentz factor and outflowing energy density (gaussian, core+power- 
law) . We find that for "on-beam" viewing geometry, with the line of sight passing through the bright 
core of the outflow, the shape of the lightcurve practically does not depend on the jet structure, 
initially following the temporal slope 2 + 8, where 8 is the spectral index. When the viewing angle is 
larger than the core, the light curve decaying slope is shallower. We discuss the implications of our 
results for the Swift data. 

Subject headings: gamma rays: bursts — gamma rays: theory 



1. INTRODUCTION 

Recent Swift observations of gamma-ray burst after- 
glows (Tagliaferri et al. 2005; Chincarini et al. 2005; 
Nousek et al. 2005; O'Brien et al. 2005) revealed very 
fast drops of X-ray flux a few hundred seconds after the 
trigger time. The rapid decay is usually also followed 
with the X-ray flares (Burrows et al. 2005; O'Brien et 
al. 2005). A possible mechanism of this phenomenon is a 
sudden decrease of jet's emissivity (see Zhang et al. 2005 
for more comprehensive discussion of possible scenarios). 
For a curved front surface of a jet the turn-off is perceived 
by an observer as a fast, yet not abrupt, decrease of flux, 
because the simultaneous drop of emissivity at different 
angles from the line of sight is perceived at different mo- 
ments. For uniform properties of the jet, and in the ab- 
sence of further complicating factors, the temporal slope 
a of a lightcurve resulting from such a "curvature effect" 
is equal to 2 + 8, where 8 is the spectral index and the 
conventions F v oc t~ a and F v cx u~ s are assumed (Ku- 
mar k Panaitescu 2000; Fan & Wei 2005; Zhang et al. 
2005; Panaitescu et al. 2005). 

In addition to the spectral properties of emitted radia- 
tion, the structure of the jet, i.e. the variation of jet prop- 
erties as a function of the angle from the jet axis, may 
influence the observed lightcurve - an issue addressed in 
many recent investigations (e.g. Zhang & Meszaros 2002; 
Rossi et al. 2002; Wei & Jin 2003; Kumar & Granot 
2003; Granot & Kumar 2003; Salmonson 2003; Zhang et 
al. 2004; Yamazaki et al. 2005). During the curvature 
effect, an observer receives radiation emitted from dif- 
ferent regions of the outflow (different 9). Therefore, an 
interesting question is whether the curvature effect de- 
pends on the unknown jet structure. The main purpose 
of this Letter is to address this question. 

2. CALCULATION METHOD 

To calculate the lightcurves we use a three dimen- 
sional code which rigorously takes into account all kine- 
matic effects that affect the observed flux (eg. Doppler 
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boost, propagation time delays) in the way described 
by Salmonson (2003). To minimize the number of pa- 
rameters involved, the curvature effect is studied for 
a short flash, with emission lasting between t- ln \ t and 
toft = I.05£i n it in the frame of the central source (we take 
?init = 10 14 cm/c, where c is the speed of light). This is 
the typical internal shock radius (Rees & Meszaros 1994), 
at which the dynamics of the outflow evolves negligibly. 
At the time t- m n, the outflow is located at the radial dis- 
tance n n it = /3((9)tinit, where /3 = v/c is the velocity of 
the outflow at the angle 9 from the jet axis, and (3 = \f3\. 
This corresponds to the choice of t = for the moment 
of ejection (see Zhang et al. 2005 for the discussion of the 
"zero time effect"), and implies that the surface of the 
structured jet is slightly non-spherical during the flash. 1 
For the emissivity in the comoving frame we assume: 



I' v {9) cx {v ohs D)- 



r 

= o 



if r(0) > to 
if r(0) < to 



(i) 

(2) 



where v h s is the observed frequency, D = r~ x (l — (3 ■ 
h)^ 1 is the Doppler factor, n is a unit vector pointing 
toward the observer, e is the energy of the outflow per 
unit solid angle, and T = (I — /3 2 )~ 1//2 is the bulk Lorentz 
factor of the outflow. The quantities e and T in eq. JIJ 
(and thereby j3 and D), all refer to the angle 9 from the 
jet axis and the redshift z — is assumed. The factor 
e/F is introduced to represent the number of electrons at 
different viewing angles, i.e. N e (6) oc e/T. 

The code is able to deal with any axially symmetric 
jet structure. After testing the code on the "top hat" jet 
case, we have made calculations for four jet structures: 



1. gaussian, with 
20? 



e(9) oc exp 



r(0) = i+(ro-i)exp 



102 
20? 



(3) 



1 Complications intr oduc ed by the ^-dependent Lorentz factor 
are discussed in Section l3.4l 
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where ro = Y(6 = 0) and 9 C is the characteristic 
width of the jet; 

2. "e-gaussian" with e(0) given by eq. but with a 
fixed r = Tq for all angles; 

3. core + power law (hereafter CPL case): 

^)*(i + ^J . rp) = i+(To-i) + pj 

(4) 

with a — 2, 6 = 2; 

4. "e-CPL" case, which differs from CPL only in that 
T is fixed (6 = in eq.®. 

Thus, we consider two qualitatively different scenarios: 
cases 1 and 3 with both T and e decreasing towards the 
edge of the jet according to (nearly) the same law, and 
cases 2 and 4 with e decreasing, but T fixed. Note that 
in cases 1 and 3 the emissivity in the comoving frame is 
uniform since e/T w constant, as is in the top-hat case. 
The variations of T influence eq. only through the 
Doppler factor D. Our calculation in the cases 1 and 3 is 
limited only to the relativistic part of the outflow (eq. |2J) ■ 

3. RESULTS 

Figs, to present our results obtained for r = 
300, 9 C = 5.7°, and 5 = 1.3. Different lightcurves, shown 
in the top panels, correspond to different viewing angles 
o bs- The solid and dotted lines correspond to the "on- 
beam" viewing (0 o b s = and 0.5# c , respectively). The 
short dashed line is for viewing along the edge of the 
bright jet core (0 o bs = 9 C ), whereas the other lines are 
for "off-beam" viewing (0 o bs > 9 C ). 

3.1. The ubiquitous slope 2 + 5 

The straight solid line in top panels marks the slope 
2 + 5 = 3.3, analytically predicted for a uniform outflow. 
It can be seen that regardless of the jet structure, ma- 
jority of modeled lightcurves follows the 2 + 5 slope for 
many orders of magnitude in flux. The only exceptions 
are the lightcurves for off-beam viewing of a jet with 
r decreasing with distance from the jet axis (gaussian 
case, Fig. and CPL case, Fig. U>). These off-beam 
lightcurves exhibit a break at a late time (t <? 30 s for the 
gaussian case, Fig. QJi, and t ^ 10 2 s for the CPL case, 
Fig.^)), after which the flux drops rapidly (on average 
following the slope ~ 3.8). The fast drop is caused by 
the lack of emission from the outer parts of the outflow 
(the jet edge given by eq. |3 starts to be visible for an 
observer) . 

The main reason for the trend to initially follow the 
2+5 slope is that the decrease of flux due to the curvature 
effect occurs on a much shorter timescale than that for 
the jet structure to take effect. For the spectral index 5 
the flux decreases by m orders of magnitude after a time 
of t crv — 10 m /( 2+s H o g , where t g is the observer time at 
which the curvature effect began. For a typical 5 ~ 1, the 
flux drops by one-order of magnitude after a short time 
i crv ~ 2t a Q. A drop of three orders of magnitude occurs 
in no more than a decade in time. On the other hand, for 



a spherical surface of the jet, 2 the observer can perceive 
the switch-off of emissivity at angle 9 measured from the 
line of sight at a time ts — [1 + (9Y) 2 ]t ff. One can see 
that the structure of the outflow must have the typical 
angular scale of the order of 3 /Y to affect the observed 
flux before 10i o ff. For r > 10 2 , the parameters of the 
outflow would have to vary strongly on a scale of one 
degree to be detectable at early time, when the flux is 
high. In both cases considered in this paper (i.e. gaussian 
and CPL), Y changes little within the core of a jet, and 
this is why a ~ 2 + 5 for 9 \,s < 9 C . The strong deviation 
from 2 + 5 starts to be pronounced only when 9 b s ~ 
9 C 2> l/To, i.e. at late times or large viewing angles. 

For a small core of angular radius 9 C = 1° the discrep- 
ancy between the calculated slope and 2 + 5 appears at 
earlier time, as can be seen in Fig. [5J 

3.2. Rebrightening 

At late observational times the lightcurves diverge from 
the 2 + 5 slope, even in the "on beam" case. An interest- 
ing feature that appears for outflows with ^-dependent 
Lorentz factor is a rebrightening, visible for the gaussian 
case (Fig. (TJi) at t > 30 s. 

The reason for this behaviour is a nonmonotonic de- 
pendence of Doppler factor D on 9. After the sud- 
den turn-off of emissivity, the observed flux decreases as 
£>( 2+<5 ) , and the Doppler factor can be approximated with 
D app = 2[Y{9)]- 1 9- 2 in the limit of 6 » 1/r. Initially, 
the Lorentz factor changes slowly (D app w 2/(Y 9 2 )), so 
that D a pp decreases with increasing 9 (and t), and the 
flux follows the slope 2 + 5. As soon as 9(t) ~ 9 C , the 
factor Y(9) starts to decrease fast which initially slows 
down the rate of the flux decrease and the discrepancy 
from the slope 2 + 5 appears. Eventually, the decrease of 
Y(9) overcomes the dependence D app cx 9~ 2 and causes 
-Dapp (as well as the observed flux) to increase with time. 

In the gaussian case (Fig.^) the flux reaches minimum 
when the observer receives radiation from 9 = 2 1 / 2 9 c , 
which happens near t ~ 30 s in Fig. da. The Lorentz 
factor at this angle is Y = 1 + (ro — l)exp(— 1) « 111 
(for T = 300 used throughout this paper). 

In the CPL case shown in Fig. the flux at t ^ 10 2 
s is nearly constant, because of the onset of the regime 
9 9 C , in which the specific structure of the outflow 
(Y(9) oc 9~ 2 ) cancels out the ^-dependence in D app cx 
T -1 ^ 2 . For the CPL case with a smaller core (Fig. EJ) 
the levelling happens at earlier time (t ~ 1 s). After 
the plateau phase the lightcurve tends to steepen again 3 
before undergoing the jet edge break. 

3.3. Role of the emission preceding the turn-off 

The lightcurves with the rebrightening effect can be 
considered as consisting of two qualitatively different 

2 The front-surface of a jet with 6-dependent Lorentz factor 
should be non-spherical because the outflow velocity is not uni- 
form. In the cases shown in Figs. Ill and l2l this effect is small and 
can be neglected (see Section 13.41 . 

3 The asymptotic slope which the lightcurve is trying to take 
on (but never succeedes because of the break) is equal to a = 
[8 + 4 — (2/n)]/2, where n > 1 is the exponent in the relation 
T(8) oc 9~ n (n = 2 in the CPL case) and the convention F v oc t~ a 
was assumed. This temporal index corresponds to the limit <C 

i/r. 
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parts: one of them has a power-law shape with a slope 
very close to 2 + 6; the other part, which includes the 
rebrightening, has a variable slope. This specific shape 
provides a good insight into the problem of how the shape 
of a lightcurve depends on the emissivity which preceded 
the turn-off. 

We have numerically tested a variety of cases, assum- 
ing different duration of emission before the turn-off, and 
various radial dependence of the emissivity. We have 
found that the flux initially drops following the famil- 
iar 2 + 5 power-law, independently on the emissivity that 
preceded the turn off. However, the shape of a lightcurve 
in the region of the variable slope (especially in the valley 
preceding the rebrightening) is very sensitive to the du- 
ration, and to the time dependence of emissivity which 
preceded the turn-off (Fig. |3J) . 

To understand this result, it is useful to consider the 
continuous emission preceding the turn-off, as consisting 
of a series of short flashes, each of which contributes its 
own lightcurve with the shape described above (power- 
law + rebrightening) but shifted horizontally with re- 
spect to each other and differently normalized. The ob- 
served lightcurve can then be considered to be a sum of 
these sub-lightcurves. In the period of time within which 
all these sub-lightcurves have the power-law shape (with 
the slope 2 + 6), the total lightcurve will have the same 
shape. In the region of the variable slope (near the re- 
brightening), the observed lightcurve takes on a different 
shape than the shape of individual components. Analo- 
gous effect happens when one calculates a photon spec- 
trum by integrating over an electron energy distribution 
having a broken power-law shape. 

Therefore, in the region of the non-power-law shape the 
lightcurves calculated for the short flash may differ no- 
ticeably from those calculated for the "sudden turn-off" 
case. Our choice of the "short flash" scenario has been 
mainly dictated by two facts: 1) The differences between 
the two cases in the non-power-law region are small if the 
emissivity preceding the turn-off increases quickly with 
time, or has a short duration. According to the internal 
shock model of GRBs (Rees & Meszaros 1994) both these 
conditions are likely to be fulfilled in reality. 2) A consid- 
erably larger number of model parameters is required to 
model the sudden switch-off case (e.g. the duration and 
the radial dependence of emissivity before the turn-off) . 

3.4. More realistic models 

For an outflow with the ^-dependent Lorentz factor, 
the standard version of a curvature effect with the turn- 
off/flash occuring simultaneously at the same radial dis- 
tance r g from the central object is hardly conceivable, 
because of the non-uniform velocity of the outflow. 

Therefore, we have also calculated lightcurves for a 
locally-short brightening of a spherical surface. The 
emission given by eq. was limited to the region 
between ri n it = Rq and 1.05i?o and started non- 
simultaneously at a ^-dependent time ti n n — r/((3{9)c) 
(which may roughly correspond to a collision of a non- 
uniform outflow with a spherical layer of material). 

Lightcurves calculated for such "non-simultaneous 
flash" at a fixed r are practically identical to those that 
we get for the simultaneous flash at the ^-dependent 
distance. This is because the additional radial delay 
A< s tr = (r/c)[l — v(9)/v(9 = 0)] associated with the 



extra curvature is much smaller than the standard de- 
lay for the spherical curvature: At sp h = (r/c)[l — cos 9]. 
In other words, the reason is the proximity of (3{9) to 1 
within the relativistic part of the outflow. 

We next proceed to discuss the case of a flash at a place 
and time that correspond to a collision of two structured 
shells: a slow one, characterized by the Lorentz factor 
profile T s {9), and a fast one [r;(#)] ejected from the cen- 
tral source St seconds later. For t = corresponding to 
the ejection of the fast shell, the collision occurs at the 
time t col (9) = 5t/{(3i{9)/l3 s {9) - 1) « 2St[T s {9)} 2 and at 
the radial distance r col (9) = Pi(9)ct co i(9) w 2c5t[T s (9)] 2 . 
Therefore, for outflows with T quickly decreasing with 9, 
as in the cases 1 and 3, the emission region is very elon- 
gated and has a pencil- like shape (r co \ oc T 2 ). Even at a 
large 9 the radiation is emitted from r co \(9) <C r co i(9 = 
0), i.e. from regions located very close to the jet axis. 

Surprisingly, however, the large additional distance 
that must be covered by photons does not produce large 
discrepancy from the 2 + 5 slope. The reason is that the 
collision at large 9 occurs at much earlier time (t co \ oc T 2 S ). 
This temporal advance [At t = i co \(0) — t co i(9); hereafter 
(0) means (9 = 0)] initially closely compensates the ex- 
cess of the spatial delay caused by the additional cur- 
vature, which is At spt = r co i(0)/c — r co i(9) cos(9)/c — 
A(0)^coi(0) — (ii{9)t co \{9) cos{9). One can see that as long 
as F; (9) ^ 1 and 9 ^> 9 C the difference between them is 
Ai spt — At t ~ At a ph, where Ai sp h = r(0)[l — cos9]/c is 
the delay for the spherical curvature. The exact value 
of the temporal advance is slightly larger than the addi- 
tional spatial delay (At t > At spt — At sp h), so that pho- 
tons emitted there are actually observed earlier than in 
the case of the spherical outflow. 

3.5. Outflow with the Lorentz factor independent of 9 

In the cases with ^-independent Lorentz factor 
(Figs.^ and^i) the lightcurves initially follow the slope 
2 + 5 and then behave in a way which is easy to in- 
terpret. For on-beam viewing (0 o bs & 9 C ) the decline of 
flux is faster than 2 + 5, because the emissivity at a late 
time (larger angles 9) is weaker. For off-beam viewing 
(#obs > 9 C ) the drop of flux becomes slower, because at 
the late time the bright core of the outflow becomes vis- 
ible for the observer. 

Unlike in the case of (9-dependent Lorentz factor, the 
problem of non-simultanous turn-off does not arise in 
cases 2 and 4. 

4. CONCLUSIONS 

We have modeled the curvature effect in structured jet 
models. For "on-beam" viewing geometry, we find that 
the temporal decay slope, 2 + 5, as is derived for isotropic 
fireballs, is ubiquitous. The deviations happen in the 
"off-beam" geometry, with the line of sight pointing to- 
wards the Gaussian or power law wing of the structured 
jet. This gives a shallower decay index than 2 + 5. 

Detailed analyses of Swift data indicate that in some 
cases, the rapid decay slopes of X-ray early afterglow 
lightcurves are somewhat shallower than 2 + 5 (O'Brien 
et al. 2005). One possible reason is that these reflect 
structured GRB jets viewed "off-beam" . 

This work was supported by a research grant at UNLV 
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Fig. 1. — Top: Lightcurves for a short flash from a structured jet illustrating the curvature effect. Columns a), b), c) and d) refer to 
the structures described as cases 1, 2, 3, and 4 in Section [3] Different lines: solid, dotted, dashed, dot-dashed, dot-dot-dot-dashed, and 
long-dashed, correspond to the viewing angles 9 y, s /9 c = 0, 0.5, 1, 1.5, 2, and 3, respectively. The straight solid line is for the analytically 
predicted slope 2 + 5. The result is for <5 = 1.3, 9 C = 5.7°, and To = 300. Bottom: Temporal slope for the lightcurves in the top panel. 
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Fig. 2. — Same as in Fig.^3 (core + power law case) but for a significantly smaller core width 6 C = 1°. 
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Fig. 3. — Comparison of a "short flash" lightcurve (solid curve) with lightcurves for the "sudden turn-off' case. The latter differ in a 
radial dependence of emissivity before the turn-off: the dotted line is for I' v ex r 3 , and the dashed one is for V v = const. Note that the 
power-law part of the lightcurve after t Q fj f» 0.02 s is not affected by the emission preceding the turn-off. 



